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Abstract. The influence of the Hall effect on the global stability of cool Kepler disks under the influence of an axial magnetic 
field is considered. For sufficiently large magnetic Reynolds numbers Rm the magnetorotational instability (MRI) exists in 
a finite interval of magnetic field amplitudes, Bmm < B < B mlix . For Kepler disks the pure MRI needs both rather high Rm 
(representing the needed electrical conductivity) as well as B mm of order 0. 1 G. The magnetic field pattern resulting from our 
\ global and linear calculations is of quadrupolar parity. For magnetic fields antiparallel to the rotation axis the Hall effect reduces 

the minimum magnetic Reynolds number by about one order of magnitude. The B mm , however, is even (sightly) increased (see 
Fig.®. 

For magnetic fields parallel to the rotation axis the Hall effect drives its own instability without the action of the Lorentz force. 
The corresponding critical magnetic Reynolds number proves to be larger with the Hall effect (Rm ~ 10) than without the Hall 
effect (Rm ~ 7) so that the Hall effect for parallel fields even disturbs the formation of MHD-instability in cool protoplanetary 
disks. If the disk is supercritical then the main result of the Hall effect for positive fields is the strong reduction of the minimum 



Or 



magnetic field amplitude which is necessary to start the instability. Observations must show whether in star-forming regions the 



rotation axis and the magnetic field orientation are correlated or are anticorrelated. If the magnetic fields are large enough then 
our model predicts the dominance of fields antiparallel to the rotation axis. 
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1 . Introduction each other close to B m i n , and they compete close to B max . Both 

boundaries of the instability range are thus reduced by the in- 

The Hall effect in protostellar disks, with their low degree of terplay of the two effects in this case, 
ionization, has recently become a subject of increasing interest 

due to its relevance to the stability and the angular momentum For ne S ative axial fields the Hal1 effect and MRI am P lif y 

. - .1 i - i rr U u u rc . , f each other when the field strength is close to B max . The Hall 

transport in the disks. The Hall effect can amplify or suppress fe mdx 

. ■ ■ . . .■ , • . u r . /u e . A/ir>T effect is destabilizing here as it transforms the value of B max 

the standard magnetorotational instability (hereafter MRI, see & °™ 

D 1U o u i inni\ a a- ,u • e t u a to higher values. The Hall effect and MRI compete, however, 
Balbus & Hawley 1991) depending on the sign of the prod- & F ' ' 

uct of angular velocity and magnetic field projections on the close to B ™<" Here ' the Hal1 effect is stabilizin S as the B ™ 

wave vector of a disturbance. The effect was found to desta- is Creased. Both boundaries of the instability range are thus 

, , i . • .. ,,,, mnn r> iu e enhanced by the interplay of the two effects in this case, 

bihze when the product is negative (Wardle 1999; Balbus & 1 F 3 

Terquem 2001; Riidiger & Shalybkov 2004). In this paper we 

shall present a study for the global stability of a differentially 

rotating disk of given (small) thickness and (low) temperature. 

As the solution of the induction equation alone, the Hall 2. Local approximation 

effect can drive its own instability. The instability does not re- 
quire rotation; it also exists for a plane shear flow. For dQ/dR < 2.1. Linearized equations 
this shear-Hall instability develops when the axial magnetic 

field is positive and vice versa. Similar to MRI, it only exists The conditions in protostellar disks were discussed by Balbus 

between a minimum field, B mm , and a maximum field, B max - & Terquem (2001). The disk material is partly ionized plasma 

For a positive magnetic field, MRI and the Hall effect amplify where ions are well linked to neutrals but electrons are not. 
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This leads to the induction equation including an additional 
term compared to the standard one-fluid MHD, i.e. 

SB 

— = Vx(« xB + H H xB - riVxB), (1) 
of 

where the second term on the right stands for the Hall electro- 
motive force with the effective velocity, «h, proportional to the 
current density, J - VxB/hq: 



Cy 



-r\ C H - 



VxB 



B 



Ve 



(2) 
(3) 



In these equations, n e is the electron number density, w ce = 
eB/cm e is the cyclotron frequency, and v e is collision frequency 
of electrons. The reason to introduce the 'Hall velocity' Q is 
that it helps in interpreting future results. Depending on the 
magnetic field direction the quantity Ch can be positive or neg- 
ative. 

The equation of motion reads 

^ + (b-V)b = --VP - V<t> + -JxB + vAu, (4) 
ot p p 

where the viscosity term is kept for numerical reasons although 
the viscosity is small for protostellar disks. We shall see that the 
stability parameters do not depend on the viscosity whenever 
the magnetic Prandtl number, 



Pm = -, 



(5) 



is below 0.1. We assume an incompressible fluid, diva = 0. 

The reference state includes a non-uniform rotation with 
the angular velocity, Q, dependent on the distance s to the ro- 
tation axis, and a uniform axial magnetic field Bo||fi. Linear 
stability or instability against small disturbances is considered. 

For a local approximation a Cartesian coordinate system 
rotating with the local angular velocity Q is used with x, y and 
z pointing in the radial, azimuthal and vertical directions (see 
Balbus & Hawley 1991; Brandenburg et al. 1995). The local 
approximation concerns perturbations whose spatial scales are 
small compared to the global scale of the disk parameters. The 
rotation law can then be approximated by the shear flow Uq = 
-e y Q.qx where q is the (constant) local shear. The linearized 
MHD equations with the Hall effect read 



dB' du' 



dB' 

xqQ, h qQ.BJe x - 

oy 



AD , n 8(VxB') 
-J]AB + T]C K = 0, 

oz 

du' du' 

— — YlLleXU - xqLl — qiluCy- 

ot ay 



\np oz p 



vAu' = 



(6) 



with dashes indicating the small disturbances. 

Considering plane waves with B',u',P' ~ exp(y? + ikz) 
leads to the dispersion relation 

((y + ijk 2 ) 2 + co H (w H - qO)) ((y + vk 2 f + 2(2- q) Q 2 ) + 



+ co\ (w 2 - 2qQ 2 + 2 (y + vk 2 ) (y + j]k 2 )) + 
+ (4-q)co 2 A oj H 0. = 0, (7) 
with the Alfven and the Hall frequencies 
w A = kB a / -v^uop, ton = rjk 2 C H - (8) 
If overstable modes are ignored, y — in gives the equation 



Pm 2 



■ ((4 - q) C H - 2qRm) + 2Ha 2 + Ha 4 = 0. (9) 



(l+C H (C H -<?Rm)) 1+2(2 -q) 



Ha 2 Rm 
Pm 



for the marginal stability separating the regions of instability 
(%(y) > 0) and stability (2v(y) < 0). Here Rm and Ha are the 
local magnetic Reynolds number and the Hartmann number 

(10) 



Rm = 



Tjk 2 ' 



Ha = 



k 2 y/rjv' 



both taken as positive-definite throughout the paper. The shear 
is obviously necessary for any instability because Eq. (|9} pro- 
vides solutions with real Ch and Ha only with finite values of 
q and Rm. 

2.2. Shear-Hall instability 



(11) 



(12) 



Consider the case where the Hall frequency l|8} is large, 
w 2 » co\ . 

Then the dispersion relation |7} reduces to 

(y + r]k 2 f + w H (con - qSl) = 0. 

The induction equation has been decoupled from the equation 
of motion, and the dispersion relation (I12> can be found from 
the induction equation alone. An instability exists if the mag- 
netic Reynolds number exceeds the minimum value 2/q. If q is 
positive it exists for positive Ch, i.e. for the external magnetic 
field parallel to the angular velocity, and vice versa (gw H must 
be positive, see J12l i). For large Rm » Rm m j n , the instability 
region is given by 

-j- < C H < ?Rm. (13) 
qRm 

The Hall parameter Q is proportional to the magnetic field. 
Similar to MRI, the shear-Hall instability exists in a limited 
range of magnetic fields. In a further similarity, the maximum 
growth rate is controlled by the local Oort A value, i.e. 



Tmax + r]k = qQ/2. 



2.3. MRI for low conductivity, Pm <sc 1 



(14) 



For the small magnetic Prandtl numbers expected for protostel- 
lar disks an appropriate scaling of the MRI parameters exists. 
Without the Hall effect (Ch = 0) Eq. @ yields the neutral sta- 
bility condition 



Rm 2 



Pm 2 (l + Ha 2 ) 2 
2 (q Pm Ha 2 - 2 + q) 



(15) 
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for MRI alone. Suppose that Pm is decreased by decreasing the 
viscosity but keeping the magnetic diffusivity finite. Then the 
Lundquist number 



1 .0 



S = VPm Ha 



Tjk 2 ' 



remains finite. At small Pm, Eq. d!5l > becomes 
S 4 



Rm z = 



2 (q S 2 - 2 + q) ' 



(16) 



(17) 



This equation does not include Pm which means that the actual 
viscosity value is not important provided that it is smaller than 
the magnetic diffusivity. We shall see in Sect. @] that a repre- 
sentation of MRI in terms of both the magnetic Reynolds and 
Lundquist numbers becomes independent of Pm when the latter 
drops below (say) 0. 1 . In a global model, therefore, a small but 
finite viscosity can be kept for numerical stability to produce 
the results which remain valid for arbitrarily small Pm when 
represented in terms of Rm and S. Hereafter, we always use the 
Lundquist number d!6i rather than Ha. 



0.6 
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Fig. 1. The growth rates, y* — y + r\k 2 in units of A-qQ./2 as 
functions of the normalized field strength for MRI alone (full 
line) and MRI modified by the Hall effect of positive fields (J3 = 
1, dashed line) and negative fields {B — -1, dashed-dotted). 
Pm = 1. 



2.4. MRI plus Hall effect for Pm = 1 

Consider the interplay of MRI and the Hall effect for the sim- 
plifying case of Pm=l. Then the marginal stability equation 
becomes 



S 2 + 2C H Rm + 1 



-Rm (2Rm + C H ) 
+ (2Rm-C H )f^-! 



1/2 



(18) 



The plus and minus signs on the RHS of this equation define 
two boundaries of the instability region for a given Rm pro- 
vided that the magnetic Reynolds number exceeds the mini- 
mum value of Rm m ; n = 2/q. 

Note that the Hall parameter and the Lundquist number 
j!6l > are both linear in the external field Bq. Their ratio 



(19) 



therefore, characterizes mainly the material. It may be taken to 
be positive for magnetic fields parallel to the rotation axis and 
negative for magnetic fields antiparallel to the rotation axis. 

For both cases the results are quite different. Consider the 
instability region for very large Rm » max(Rm m j n ,y8 _1 J. For 
positive 8 it is 



1 



S < Rm((2 9 +^ 2 )' /2 -jSj. 



sqRm --W- -> -r (20) 

The lower limit is the same as in Jl 31 . B m i n for positive B is 
thus controlled by the shear-Hall instability. The lower bound 
is small compared to S m m = yj(2 - q)/q expected from Eq. 
fP71 for MRI. The Hall effect for positive B Q amplifies MRI for 
the fields close to Z? m j n . The effect is, however, stabilizing close 
to B max because the upper limit in ^20^ decreases with B. 
The instability domain for negative B is 



|(2-4)Rm<S < ^2q + B 2 + \B\jRm. 



(21) 



Both bounds increase with \B\. Here the Hall effect is thus desta- 
bilizing close to B max (Wardle 1999; Balbus & Terquem 2001) 
and stabilizing close to B m i n . 

In Fig. [2 the interaction of MRI with the Hall effect is 
demonstrated by the growth rates derived from the dispersion 
relation @ for Pm = 1 . The solid line gives the typical growth 
rate profile for the MRI alone. For positive fields the profile is 
moved by the Hall effect to the left and for negative fields the 
profile is moved to the right. At the strong-field limits only the 
negative Hall effect is thus destabilizing. In opposition to that 
the positive Hall effect even disturbes the instability and may 
not help to produce the desired turbulence in cool disks. 

3. A global model 

The above results for growth or decay of small disturbances 
taken from the local dispersion relation have been known since 
the papers of Wardle (1999) and Balbus & Terquem (2001). 
Contrary to the local considerations we shall present in the fol- 
lowing the results for global Kepler flow models with Ohmic 
dissipation and the Hall effect included. The progress of such 
an approach is to find exact values for the critical Reynolds 
numbers, the minimal and maximal magnetic amplitudes and 
the global geometry of the resulting instabilty pattern. As the 
limiting magnetic fields prove to be surprisingly high (and 
rather sensitive to the Hall effect) information about the min- 
imum magnetic fields is important for the discussion of the 
(magneto-)hydrodynamics of cool protostellar disks. 

Our global model differs from that of Kitchatinov & Mazur 
(1997) only by the inclusion of the Hall effect. The model con- 
cerns a rotating disk of constant thickness, 2H, threaded by a 
uniform axial magnetic field. The rotation axis is normal to the 
disk, and the angular velocity, £2, depends on the distance, R, 
to the axis. This dependence is parameterized by 



Sl(R) = n Cl(R) 



(22) 
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with 



1+1* 



3«/2 \ 



-1/n 



(23) 



This profile describes almost uniform rotation at small dis- 
tances R <sc Rq, which smoothly transforms to the Keplerian 
law, Q ^ Q.o(Ro/R) i/2 , for large distances R » We use 
n — 2 in Eq. d23i and Ro/H = 5 for the aspect ratio and 
div u — as above. The pressure is excluded by curling Eq. 
10} . This yields for the vorticity, to — V x u, the relation 



3°* „ , 

— = V x (u x to 

dt 



J x B/p) + vAco. 



(24) 



The Eqs. Q and J24> are about the rotation 1221 and the uni- 
form axial field, Bo = fio<^ (2* is me un it vector along the 
rotation axis), and the normalized variables 



(25) 



b = B'/Bo, j = ^-J', u = u'KHClo), to = oZ/Qo, 

DO 

for the disturbances are introduced. This leads to four basic di- 
mensionless parameters among which the Hall parameter, Ch, 
and magnetic Prandtl number, Pm, are defined by (0 and 
but magnetic Reynolds number and the Lundquist number now 
are 

CloH 2 „ B ( )H 

(26) 



Rm 



/WW 



These parameters control the equation system for the normal- 
ized disturbances 



Rm V x (R&(R) c x b - e z x u) - C H (e z ■ V) }" + Ab, 



db 

dt 

do) I - k 2 

— = Rm V x RQ.(R) e^xto : e xu 



(27) 



(28) 



+ (e, • V) j + PmA<£>, 

Rm 

where k is the normalized epicycle frequency 

2 _ 2 Q d(^ 2 Q) 

K ~ R dR ' 

and time and distances are normalized to the diffusion time, 
H 2 /tj, and the disk half-thickness, H. Here only the stability of 
the rotation law is considered, the accretion flow which is con- 
nected to this nonuniform rotation via the viscosity is neglected 
(see Kersale et al. 2004). 

The boundary conditions on the disk surfaces are (i) stress- 
free for the flow and (ii) pseudo-vacuum conditions for the 
magnetic field fluctuations, i.e. 



e, x b = 0. 



(29) 



The solutions are required to be regular on the rotation axis and 
to vanish at infinity. 

The linear stability analysis with b, to ~ exp(yf) leads to an 
eigenvalue problem for the equation set J27l > which has been 
solved numerically. A new variable, y, 



R/Rq 
1 +R/R ' 



0<y< 1, 



(30) 



UNSTABLE 



Pm = 0.001 
Pm = 0.01 
Pm = 0.1 
Pm=1 



1 

s 



Fig. 2. Neutral stability lines of the MRI for various magnetic 
Prandtl numbers (Ch = 0). For Pm < 0.1 the lines are almost 
independent of Pm. 

has been introduced transforming the infinite disk to a finite 
domain. Rq is the turnover radius in Eq. (I23i . A uniform grid 
in y was applied which corresponds to a non-uniform grid in R. 
The eigenvalues and eigenvectors are computed by the inverse 
iteration method. 

The system ( I27l > allows two types of solutions with differ- 
ent symmetries about the disk midplane. One of the symmetry 
types combines a symmetric magnetic field with an antisym- 
metric flow field. The notation, Sm, will be used for this type 
of eigenmodes, where 'm' is the azimuthal wave number, i.e. 
SO represents an axisymmetric mode, SI defines the nonax- 
isymmetric mode with m = 1, and so on. The other symmetry 
type combines antisymmetric magnetic field with symmetric 
flow. The notation Am is used for the eigenmodes of this type 
of symmetry. 

The primary goal of the linear theory is to define the stabil- 
ity boundary in parameter space which separates the region of 
stable perturbations with negative or zero real part of y from the 
instability region with (exponentially) growing perturbations. 
The stability map strongly depends on the symmetry type of 
the excitation. 

4. Results 

4.1. MRI for low conductivity (small Pm) 

Consider first with Ch = the MRI alone. Figure [2] shows 
the neutral stability lines for Prandtl numbers decreasing from 
Pm=l to smaller values. The stability is almost independent 
of the magnetic Prandtl number for sufficiently small Pm, say 
Pm < 0. 1 . The MRI characteristics computed with moderately 
small Pm remain valid for arbitrary small Pm. This finding may 
be important for predicting MRI parameters for laboratory ex- 
periments and protostellar disks where magnetic Prandtl num- 
bers are very small. Hereafter, we fix Pm = 0.01. 

The minimum value for Rm for instability in the small Pm 
regime is 

(31) 
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Fig. 3. Stability diagram for the shear-Hall instability (positive 
magnetic field, i.e. parallel to the rotation axis; only the induc- 
tion equation is solved). Lines are labeled by their symmetry 
types. The axisymmetric quadrupolar mode (SO) is preferred. 



For Reynolds numbers large compared to this value the insta- 
bility exists within the magnetic range 

1 < S < 0.5 ■ Rm. (32) 

Cool protostellar disks may not reach such values of the mag- 
netic Reynolds number. On the other hand, for protostellar 
disks the minimum condition for the magnetic field (S 1) 
proves to be very stringent (in great contrast to the MRI in 
galaxies, see Kitchatinov & Riidiger 2004). With the values 
given below for 77 (10 15 cm 2 /s), density (10~ 10 g/cm 3 ) and for 
a disk height of 0. 1 AU a minimum field of almost 0.1 G is 
needed to fulfill S = 1 . It is thus tempting for different reasons 
to probe the shear-Hall instability for protostellar disks which, 
however, only exists for magnetic fields parallel to the rotation 
axis. 

4.2. Shear-Hall instability 

Another extreme exists for magnetic fields parallel to the rota- 
tion axis. The Hall effect for one sign of the magnetic field in 
connection with differential rotation can form its own instabil- 
ity as shown by the solution of the induction equation alone. In 
the present model this shear-Hall instability can be found for 
small Lundquist number, S <c Ch, i.e. large /3 of Eq. dl9> . In 
this case the Lorentz force in (0} can be neglected so that MRI 
is excluded. 

Figure [3] shows for (3 — > 00 the stability map for modes of 
three basic symmetry types, ie. SO, SI, and AO. The quadrupo- 
lar axisymmetric modes SO are preferentially excited. The neu- 
tral stability lines for the other symmetry types lie completely 
inside the unstable region of SO. It is one of the advantages 
of our approach that the equatorial symmetry results from the 
computations rather than needing to be prescribed. Salmeron & 
Wardle (2003) in their extensive analysis of the vertical struc- 
ture of accretion disks used those symmetry conditions at the 
equator valid only for AO-fields. When the system becomes 




Fig. 4. Magnetic field vectors in the disk midplane for 
marginally stable SO-modes at Rm = 100 for the minimum 
(left) and maximum (right) Lundquist numbers, resp. The bro- 
ken circle shows the turnover radius sq of the rotation law. 



unstable against AO-disturbances it is thus already unstable 
against SO-perturbations. 

The minimum Reynolds number for all the modes is 

Rm^n - 10. (33) 

For large Rm, the left branch of the neutral stability line for 
SO-modes approaches the relation 

C H - 5/Rm, (34) 
while the right branch is close to 

C H - 0.2 • Rm. (35) 

Note that for large magnetic Reynolds number Rm the mini- 
mum Ch necessary for instability becomes infinitely small. 

The maximum fields allowing for the axisymmetric SO- 
modes and the nonaxisymmetric Sl-modes in Fig.|3]are almost 
equal. This property might be important for the dynamo theory 
with respect to the Cowling theorem. 

Equations (1331 - d35l correspond to dl 31 obtained with the 
local analysis. Quantitative differences are partly due to the dif- 
ference in definitions of Rm between local and global calcula- 
tions. The differences are actually quite small if the transfor- 
mation rule k — > n/H between local and global formulations is 
used. 

Figure 0] shows that the pitch angle for marginally stable 
disturbances is indeed very small for the minimum field pro- 
ducing the instability. The angle increases to about n/2 for the 
maximum field. This tendency for the pitch angle can also ex- 
plain the results for nonaxisymmetric disturbances. For pitch 
angle close to n/2 the azimuthal structure is not significant. 
Accordingly, the maximum fields for the S 1 and SO modes are 
roughly the same (Fig. Small pitch angle, however, would 
mean a tight winding with a small radial scale for SI modes 
whose instability is then precluded by diffusion. This is why the 
instability region of SI modes in Fig. [3] exhibits such a sharp 
and almost vertical boundary on the weak-field side. 

The orientation of the magnetic field in relation to the ro- 
tation axis plays an important role in the interplay of differen- 
tial rotation and the Hall effect. Now the two cases of paral- 
lel and antiparallel magnetic fields are considered solving the 
complete equations. We find that the MRI of the Kepler flow is 
very differently modified by the Hall effect. 
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Fig. 5. The same as in Fig.[3]but for the full equation system. 
The Hall effect increases the minimum Rm. The solid curve for 
ft — > oo is identical to the lower curve in Fig. [3] For increasing 
Rm the lower limit for the magnetic field becomes smaller and 
smaller. 
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Fig. 6. Stability diagram for MRI modified by the Hall effect 
for magnetic fields antiparallel to the angular velocity. Note 
that the Hall parameter ft for real cool protostellar disks is of or- 
der unity. The Hall effect shifts the instability interval to larger 
Lundquist numbers, i.e. to higher magnetic fields. 



4.3. Positive (parallel) fields 

For positive /3 the solution is located between the two realiza- 
tions given in Figs.|5]and|3] The minimum magnetic Reynolds 
number moves from the value 7 for MRI to about 10 for the 
shear-Hall instability. The Hall effect does thus not support the 
instability of the cool Kepler flow; for all positive ft the mini- 
mum remains between 7 and 10 (Fig. |5j. 

Figure|5]also shows that the minimum of the possible mag- 
netic fields is fixed by the Hall effect. We find Rm ~ 1 1 (ft ■ S) 
so that S ~ I lift - Rm) which for large Rm is smaller by orders 
of magnitude than S ~ 1 taken from Fig.|5]for negative ft. 

Our main interest also concerns the results for low electri- 
cal conductivity, i.e. for small Rm. The minimum Rm moves in 
Fig-Elfrom 7.1 to 10, i.e. it increases opposite to the desired 
trend. If the results for parallel and antiparallel magnetic fields 
are compared, one of the differences is the opposite trend for 
the minimum magnetic Reynolds number which only for neg- 
ative fields (antiparallel to the rotation axis) is reduced by the 
Hall effect. It is also important that for large enough Rm the 
minimum magnetic fields for the instability strongly differ for 
both the magnetic orientations. 

4.4. Negative (antiparallel) fields 

Figure[6]shows the stability diagram for negative Bo f° r which 
the shear-Hall instability does not exist. The Hall effect, how- 
ever, increases both instability limits of the field amplitudes. 
The Hall effect thus destabilizes for strong fields close to B max 
and stabilizes close to B m ; n . 

On the other hand, the absolute minimum of the magnetic 
Reynolds number Rm is reduced. The necessary electrical con- 
ductivity of the gas is (slightly) reduced by the Hall effect. 
From numerical arguements, we can present only the results 
for a small Hall parameter ft which is still too small by one or- 
der of magnitude (see Eg. 1401 below). A massive reduction of 
the critical magnetic Reynolds number (i.e. the necessary elec- 



Table 1. Ratio of the global magnetic to kinetic energy at min- 
imum magnetic Reynolds number for a sequence of the Hall 
parameter ft. 



p 


-0.3 


-0.1 





0.1 


0.33 1 


^rnag / ^kin 


0.77 


1.20 


1.52 


1.94 


3.32 11.4 



trical conductivity) is thus expected from the Hall effect for 
antiparallel magnetic fields. It seems indeed to be possible to 
use the Hall effect to realize the MRI also for the rather low 
electrical conductivity of cool disks. 

4.5. Energy relation and angular momentum transport 

Linear computations cannot provide energy values. It is pos- 
sible, however, to compare their global magnetic and kinetic 
energies. These ratios of the total (volume integrated) energies 
are given in Table^for a sequence of values of the Hall param- 
eter ft. The total energy is normally dominated by its magnetic 
part in agreement with other simulations of MRI (cf. Stone & 
Norman 1994; Brandenburg et al. 1995; Sano & Stone 2002). 
The energy ratio which we obtain is, however, decreasing and 
even drops below unity when the Hall parameter becomes more 
negative. A similar trend can also be found in the simulation 
results of Sano & Stone (2002). The angular momentum al- 
ways flows outwards, and it is dominated by the Maxwell stress 
(Balbus &Hawley 1991; Brandenburg et al. 1996). The relative 
contribution of the Reynolds stress also increases with decreas- 
ing ft, in accordance with the results of Sano & Stone (2002, 
their Table 4). 

5. Discussion 

Ionization of protostellar disks material should be extremely 
low (Gammie 1996). The electrical conductivity is then con- 
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trolled by electron-neutral collisions and the magnetic diffusiv- 
ity 



/; = 234 — T 1/2 cirr s 



(36) 



is inversely proportionate to the ionization ratio (cf. Balbus & 
Terquem 200 1 ) where « e and n are the number densities of elec- 
trons and neutrals. Then the magnetic Reynolds number J26> 
can be estimated with Eq. ( I36l > as 

Rm = 2aO>^)(^n^r 1 UU\ (37) 
In/UOOK/ \0.1yr/ \0.1AU/ 

where r rot is the rotation period at the distance so = 5H. 

The ionization fraction, « e /n, is very uncertain. We estimate 
first its value required for instability. Without the Hall effect our 
results suggest Rm 10 for instability. Assuming numerical 
values in Eq. fl37i as representative, one finds 



— a 10 

n 



■14 



(38) 



as the critical value. Even such a low ionization is problem- 
atic for protostellar disks (Stepinski 1992; Gammie 1996). 
Collisional ionization is inefficient for T < 10 3 K. The cos- 
mic ray, however, can provide ionization fractions of n e /n 
10~ 12 if the column density S does not exceed 10 2 g cirT 2 
(Umebayashi & Nakano 1981; Gammie 1996). With this value, 
Eq. O for thick disks (H a 0.1 AU) provides Rm^lO 3 but 
this value reduces to 10 for thinner disks (H a 0.01 AU). 
The Hall parameter Q can be estimated as 



/ T y l/2 ( n y l B 

CH " 20 llooKj lio^wi Tg 

and B reads B of Eq. iV)\ , reads 



H V 



1AU 



(39) 



(40) 



With a particle density of n a 10 14 cm 3 , this relation leads to 
the value of 



a 2, 



(41) 



so that the Hall effect should indeed be very important. With 
such large values one can take from Fig.[3]that the pure shear- 
Hall instability dominates for positive fields. In Fig. [7]the nu- 
merical results are summarized for the dependence of the crit- 
ical magnetic Reynolds number on the Hall parameter B. For 
negative B of order unity it is obvious that the critical magnetic 
Reynolds number is reduced by at least one order of magnitude 
by the Hall effect. 

The boundaries of the instability domain also strongly de- 
pend on the field orientation. Only for antiparallel fields does an 
instability exist for Rm < 7. If the magnetic Reynolds number 
exceeds O(10) then both magnetic orientations lead to instabil- 
ity but for rather different magnetic amplitudes. It is 



0.1 G < B < 10 G 

for antiparallel fields and 

0.001 G < Bo < 1 G 



(42) 



(43) 




■1.0 



-0.5 



0.0 



0.5 



1.0 



Fig. 7. The computed minimum magnetic Reynolds number as 
a function of the Hall parameter. B = denotes MRI. The re- 
sults for B < -0.3 are not yet known. 

for parallel fields, both taken for Rm = 100. Note that the mag- 
netic fields of meteorites vary between 1 G and 10 G, close to 
the upper limits of the above equations. This coincidence sug- 
gests that the instability may drive a dynamo which saturates 
when the field is amplified to the upper boundary of the insta- 
bility interval. Also, the minimum magnetic fields allowing the 
instability prove to be rather strong. If the external magnetic 
fields never exceed 0.1 G in amplitude then only the Hall ef- 
fect of parallel fields would lead to the instability necessary to 
remove the angular momentum in the Kepler disk but then the 
ionization must be high enough. If it is not, then the Hall ef- 
fect for antiparallel fields is needed but the magnetic field must 
exceed 0.1 G in this case. It this clear from such considera- 
tions that in dense and cold globules of molecular clouds there 
are severe limitations to the necessary transport of angular mo- 
mentum by magnetic instabilities. 
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